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ABSTRACT 

We predict the Tully-Fisher (TF) and surface-brightness-magnitude relation for disk galaxies 
at z — 3 and discuss the origin of these scaling relations and their scatter. We show that 
the variation of the TF relation with redshift can be a potentially powerful discriminator of 
galaxy-formation models. In particular, the TF relation at high redshift might be used to break 
parameter degeneracies among galactosynthesis models at z = 0, as well as to constrain the 
redshift distribution of collapsing dark-matter halos, the star-formation history and baryon 
fraction in the disk, and the distribution of halo spins. 
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1 INTRODUCTION 

With the advent of numerous lOm-class telescopes, adaptive op- 
tics, and high-z galaxy surveys, it is becoming possible to learn 
about the detailed properties of galaxies at high z. In the future, we 
may expect improved IR sensitivity with NGST and 30m or even 
100m telescopes, and this should usher in an era of high-precision 
kinematic and photometric studies of high-z galaxies. Disk galax- 
ies are particularly important to our understanding of galaxy for- 
mation and evolution, as they are believed to undergo a relatively 
smooth formation process, and possibly serve as building blocks in 
the formation of other galactic systems through mergers. 

Many authors have investigated galactosynthesis models (e.g., 
Dalcanton, Spergel & Summers 1997; Mo, Mao & White 1997; 
limenez et al. 1998; Somerville & Primack 1999; van den Bosch 
2000; Avila-Reese & Firmani 2000; Firmani & Avila-Reese 2000; 
Navarro & Steinmetz 2000) in which the properties of disk galaxies 
are determined primarily by the mass, size, and angular momenta 
of the halos in which they form, and which may contain the ef- 
fects of supernova feedback, adiabatic disk contraction, cooling, 
merging, and a variety of star-formation (SF) recipes. In an earlier 
paper, we (Buchalter, limenez & Kamionkowski 2000; hereafter 
BIK), investigated a variety of galactosynthesis models with real- 
istic stellar populations and made multi-wavelength predictions for 
the Tully-Fisher (TF) relation. With reasonable values for various 
cosmological parameters, spin (A) distributions, formation-redshift 
(zf) distributions, and no supernova feedback, we could produce an 
excellent fit to the local TF relation at all investigated wavelengths 
(B, R, I, and K), as well as to B-band TF data at z — 1, and 
to the surface-brightness-magnitude (fi-M) relation locally and at 
z = 0.4. These successes suggest that our simplified approach cap- 
tures the essential phenomenology, even if it leaves out some of the 
details of more sophisticated models. 

In this paper, we investigate the high-z TF predictions for our 



model. We examine the key factors impacting the TF relation at 
high z, and demonstrate that degeneracies among parameters in 
galactosynthesis models at z = can be disentangled at high z. 
Specifically, information about the Zf and A distributions, as well 
as the SF history and baryon fraction in the disk, can be gleaned 
from the high-z TF relation. 



2 THE MODELS 

Here, we briefly review the main ingredients of our model and refer 
the interested reader to BIK for a comprehensive description. We 
use the spherical-collapse model for halos, the distribution of halo- 
formation times from the merger-tree formalism (Lacey & Cole 
1993; Lacey & Cole 1994), and a joint distribution in A and v, the 
peak height (Heavens & Peacock 1988; Catelan & Theuns 1996). 
Halos are treated as isothermal spheres with a fixed baryon frac- 
tion and specific angular momentum, and their embedded gaseous 
disks, assumed to form at virialization, have an exponential den- 
sity profile. We implicitly assume that the halos of spiral galaxies 
form smoothly, rather than from major mergers (BIK; Eisenstein & 
Loeb (1996)). An empirical Schmidt law relates the star-formation 
rate (SFR) to the disk surface density (Kennicutt 1998). A Salpeter 
initial mass function, a prescription for chemical evolution, and a 
synthetic stellar-population code (limenez et al. 1998) provide the 
photometric properties of disks at any z. 

The model is defined by cosmological parameters and by the 
time when the most massive progenitor contains a fraction / of the 
present-day mass, when a halo is defined to form. We found that 
excellent agreement with current data was obtained by our 'Model 
A', a COBE-normalized ACDM cosmogony (T2 = 0.3, h = 0.65, 
Q b h 2 = 0.019) with / = 0.5. The TF relation in this model re- 
lied on both halo properties and SF in the disk. The local TF scat- 
ter arose primarily from the Zf distribution, and secondarily from 
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chemical evolution and the v-\ anticorrelation. In this model, disk 
formation occurs primarily at 0.5 ^ z ^ 2 and the TF slope steep- 
ens and the zero points get fainter from z = to z — 1. Moreover, 
the amount of gas expelled from or poured into a disk galaxy in 
this model is relatively small and the disk and halo specific angular 
momenta are equal. 

A suite of other models that give good fits to the observations 
at low z, can also be found. To illustrate, we examine two alterna- 
tives that yield comparably good fits to the slope, zero-point, and 
scatter of the TF relation at z — in B, R, I, and K. The first 
is a CDM model with flo = 1, h = 0.65, constant metallicity 
given by the solar value, A = 0.05 for all disks, a power-spectrum 
amplitude as = 0.5, an empirical shape parameter V = 0.2, and 
/ = 0.5. High-^o models generally produce disk galaxies too faint 
to lie on the local TF relation. To compensate for this, we assume 
fif,/i 2 = 0.045, and thus term this the 'high-fib' model. 

Our second alternative is a ACDM cosmogony like Model A, 
but with metallicity held constant at the solar value, A = 0.05 for 
all disks, and / = 0.9. This results in a narrow distribution of 
formation times peaked around z ~ 0.2 (for L*-type disks), with 
appreciable ongoing formation today and almost no halos form- 
ing earlier than z = 1 (cf., Figure 1 of BJK). We thus denote this 
the 'low-2/' model. This model will produce extremely young and 
bright disks (cf., Figure 5 of BJKjQ To compensate for this, we re- 
duce the efficiency of SF in our Schmidt-law prescription by 50% 
[see equations (2)-(4) in BJK]. This lower SF efficiency is plau- 
sible given a possible inconsistency between the Model-A predic- 
tions for the absolute disk gas fractions and those observed (see 
BJK for details). 

Figure |j| depicts the z = predictions for the A, high-f^, 
and low-2/ models, in the B, R, I, and K bands. The open sym- 
bols represent extinction-corrected data from Tully et al. (1998). In 
each plot, the data are fit to the corresponding model, and p gives 
the probability for the data given the model. Since the data have ex- 
cluded spirals that show evidence of merger activity or disruption, 
we exclude from our predictions those galaxies with B — R < 0.3. 

Each of the three models yields a reasonable fit to the slope 
and normalization of the TF relation in all wavebands. Moreover, 
the predicted scatter in the B, R, I, and K bands roughly agrees 
with the observed values of 0.50, 0.41, 0.40, and 0.41, respectively, 
with the largest discrepancies occurring in the bluer bands, for rea- 
sons discussed in BJK. 



3 THE TF RELATION AT HIGH Z 
3.1 TF Scatter for a Fixed Mass 

Figure ^ shows predictions for our A and high-fii, models at z — 3. 
At this z, the low-2/ model fails to form disks. In the left-hand 
panels, the points for a given mass scatter in a nearly vertical direc- 
tion in contrast to the more horizontal scatter in the z = results 
(see BJK Figure 13), particularly in redder bands. For z = 0, the 
TF scatter arises primarily from the spread in Zf, with the earliest- 
forming objects of a fixed mass having higher circular velocities. 
This is because for an evolved galaxy that has converted most of its 
gas to stars, the light from giants serves as a faithful tracer of the 
total mass, regardless of the precise age of the galaxy. As one goes 
either to higher masses in a given waveband, or to bluer bands at a 
fixed mass, the scatter in the local TF relation for a fixed mass does 

* We note that bulges, however, may be present at high z. 



begin to curve upward (i.e., later-forming disks appear brighter) to 
form a 'magnitude peak,' since in both cases one becomes increas- 
ingly biased towards younger, brighter systems. 

At * = 3, the scatter still owes mainly to the spread in Zf, 
and Figure ^ still indicates that at 2 = 3, the very earliest-forming 
objects do possess somewhat higher circular velocities for a given 
mass. (Note in the left panels the 'peak' of points formed at the 
high- 14 end of the four scatter plots for the different masses.) The 
galaxies that populate the TF relation at z J> 3, however, are the 
small minority that constitute the exponentially decaying high-2 
tail of the adopted Zf distribution (BJK Figure 1). Thus, the forma- 
tion of these halos spans a very short epoch, so that halos of a given 
mass will have roughly the same 14, resulting in the tight horizontal 
scatter for each mass. The large vertical scatter arises from the fact 
that, though these galaxies will eventually form the lower envelope 
of the present-day TF relation, comprising the oldest and (in bluer 
bands) faintest objects of a given mass, at z = 3 even the least 
massive of these galaxies are typically no older than a few hundred 
million years, and the more massive systems younger still. Thus, 
systems on the 'young' side of the M peak have formed so close 
to the epoch of observation that they simply have not had sufficient 
time to build up appreciable stellar populations; given their high 
gas fractions and SFRs, even a small increase in their age results in 
a substantial increase in their stellar luminosity. 

At high z, disks with circular velocities 200-500 km s _1 
(masses ~ 10 13 Mq) correspond to rare fluctuations. Still, at these 
z, one does expect to find several such objects within the horizon 
in a ACDM cosmogony, and since these will be the brightest ob- 
jects in the sample, we include them in the discussion here. We also 
reiterate that the disks that form at z i; 3 are not the dominant pop- 
ulation of disks at z = 0, as most disks at z = probably formed 
at 0.5 <J 2 <J 2 (see, e.g., BJK). Disks forming at higher z will de- 
lineate the lower envelope of the TF relation at 2 = 0; for a given 
mass, these older, higher- Vc systems will be fainter since they will 
be the first to run out of fuel. 

3.1.1 Overall TF Relation 

The trend of steeper TF slopes and fainter zero-points in going from 
2 = to z — 1 is seen to continue here out to 2 = 3. Specifically, in 
the B band, the 2 = 3 TF relation goes from ~ 2 mag brighter than 
its local counterpart at log 214 ~ 3 to less than 1 mag brighter at 
log 214 = 1.7. In K, the high-2 prediction is about 1 mag fainter at 
log 2V C = 1.7 and comparable to the local prediction at log 2V C = 
3.0. This is because for fixed 14, halos at high z will be younger, 
smaller, and have lower A than halos today. In massive systems 
most of the SF activity takes place at early times. Thus, younger 
systems will be significantly brighter, particularly in bluer bands, 
even though much of their gas remains unprocessed. Smaller, low- 
14 systems, on the other hand, do not convert gas to stars as rapidly 
and are thus fainter than their local counterparts. 

The overall A-model B-band scatter in the high-2 TF relation 
is comparable to that at z = 0. In K, however, the scatter at 2 = 3 
has increased by a factor of two. Moreover, the variation in the TF 
scatter among different wavebands is not as great at high z as was 
seen in the local case. This is because (a) disks at z > 3 are likely 
to be in a starburst-dominated phase since there is very little cosmo- 
logical time for even single stellar populations to decay; (b) at such 
early times the spread in mass of gas turned into stars is very large 
for different Zf and A [see eq. (7) of Heavens & Jimenez (1999)]. 
Therefore, it seems natural that all bands reflect the maximum pos- 
sible spread. By contrast, there is little change in TF scatter from 
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Figure 1. Predicted B-, R-, I-, and K-band TF relations at z = for the A model, high-Si;, model, and low-2 1 model. The scattered dots show the results for 
fixed masses of 10 10 , 10 11 , 10 12 , and 10 13 Mq. The solid lines are the best-fitting TF relations, with zero points and slopes given by a and b, respectively, 
while the dashed lines show the 1<t scatter, denoted in each plot by <j. 



z = to z = 1, as the majority of galaxies are sufficiently evolved 
by z = 1 (Figures 13 and 14 in BJK). Furthermore, the scatter ob- 
tained here should be considered as an upper limit since, in practice, 
M-limited surveys may not detect arbitrarily young, and therefore 
faint, objects, particularly at lower masses^] However, different re- 
alizations of a given model can generally produce results that vary 
by < 0.1 mag in the TF scatter. Given the simplifying assumptions 
of the model, attention should be paid to the overall magnitudes 
of and relative differences between the slopes, normalizations, and 
scatters, rather than their precise values. 



3.1.2 Spin-Peak-Height ( Anticorrelation 

In BJK, we also noted the slight reduction in TF scatter at z — 
in Model A arising from the existence of a weak anticorrelation 
of A with v, and thus with Zf. This reduction in TF scatter was 
only about 0.15 mag in B and 0.05 mag in K for Model A, and 
lesser or no effect was seen in some other possibly viable models. 
To investigate the effect of the joint \-v distribution, we re-ran the 

t Our models employ a minimum time step of 10 7 yr — corresponding to 
the timescale for giant-molecular cloud destruction (e.g., Jimenez et al. 
2000) — that effectively dictates the faintest objects at a given mass, and 
therefore the maximum possible scatter. 



models in Figure g with A = 0.05 for all halos. This produces 
TF predictions with shallower slopes and brighter zero-points. This 
is due to the fact that the A = 0.05 value is higher (lower) than 
the full model would typically assign to high-mass (lower-mass) 
galaxies, which thus become correspondingly fainter (brighter) in 
the fixed-A model. More interestingly, using a fixed A at high z 
reduced the overall TF scatter by about 0.2 mag in all bands at z — 
3, as compared with using the joint distribution in A and Zf. This 
is in direct contrast to the z = case, where a fixed A produced 
a comparable or larger scatter, as described above. This is due to 
the vertical nature of the scatter in the TF for a fixed mass at high 
z. In this case, higher- A systems are still fainter, but their lower z; 
(already confined to only a narrow range) can only further scatter 
them to lower luminosities, not to lower V c as is in the local case, 
thereby increasing the scatter. This suggests that the scatter in the 
TF relation at high z could be used to gauge the strength of the 
actual X-Zf anticorrelation. 



3.2 The Surface-Brightness-Magnitude Relation 

Figure |3| shows the surface-brightness-magnitude relation at z — 
0.4 and z = 3 for the A and high-£7f, models, The dots are from 
2dF (Driver & Cross 2000), and the open triangles from the HDF 
(Driver 1999); the latter have z ~ 0.4. Since our model does not 
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Figure 2. TF predictions of the A model and high-Hf, model at z = 3. The low-2 1 model predicts essentially no disks to have formed at z = 3 at these scales. 



include a bulge component, the predicted /i may underestimate ob- 
served values if the disk and bulge components cannot be resolved. 

For Model A, the slope of the relation at z = 3 is slightly 
shallower than the local relation, with [i nearly identical to the low- 
z case for Mb — 51og/t > —16 and fainter by about 1 mag for 
Mb — 5 log ft = —24. Thus, at these scales, the trend towards 
higher luminosities and smaller disk scale lengths is offset by cos- 
mological fi dimming. A more drastic contrast is seen in the high- 
fli, model, where the predicted /i at the bright end varies by al- 
most 4 magnitudes. This owes largely to substantial brightening 
that occurs at late times as the high baryon fraction is converted 
into starlight. 

For both models, we find that the scatter at z > 3 is much 
tighter than for the local case because the spread in the radii of disks 
at z > 3 is smaller than at low z. For example, at z ~ 0, disk sizes 
span a range roughly between 1 and 100 kpc, while at z ~ 3 they 
range only from about 0.5 to 15 kpc in size. On the other hand, the 
respective spreads in M differ by no more than a factor of two. The 
small spread in the relation then reflects the small spread in disk 
sizes predicted in our adopted structure-formation scenario. Thus, 
it will be interesting to learn whether this narrow spread is in fact 
observed at high z. If a much wider spread is observed, this might 
indicate that either disks at z > 3 are more substantially affected 
by merger activity or that they are simply larger; i.e., the spread in 
sizes predicted by our simple isothermal spherical-collapse model 
with angular-momentum conservation is simply too narrow. Note 
that the fi predictions in general are more sensitive to assumptions 



regarding A than are TF predictions, thus providing a strong com- 
plementary constraint on this aspect of the model. 

3.3 Addressing Possible Degeneracies 

We now address how high- 2 predictions can be used to distinguish 
among galactosynthesis models that are degenerate at z = 0. The 
right panels of Figure [| show TF predictions for the high-fib model. 
At these mass scales, the low-2/ model predicts essentially no disks 
(and thus no TF relation) to be present at * = 3, certainly a drasti- 
cally different prediction than that of the other models, despite their 
close agreement at z = 0. 

The precise slopes and zero-points for the A and high-fib 
models differ more appreciably than in the z — case, with the lat- 
ter model yielding predictions that are significantly fainter at these 
scales. This owes both to the different normalizations of the fluc- 
tuation spectrum in the two models, and to the fact that the lower 
SF efficiency in the high-fif, model results in fainter disks at early 
times. More importantly, the predicted TF scatter in the high-fi;, 
model is substantially less than that of the A model. This is be- 
cause (a) a higher value of flo implies a narrower age range for 
disks forming at z > 3; and (b) the scatter due to a spread in A is 
excluded in this model. 

Alternatively, one could envision a high-z/ model in which 
disks were somehow brightened at later times, perhaps through late 
infall of fresh gas. In this case, there might be galaxies that are rea- 
sonably evolved, even at these high z, so that many galaxies would 
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Figure 3. Predicted surface-brightness-magnitude relation for disk galaxies 
at z = 0.4 and z = 3 for the A model and high-fi;, model. The solid lines 
are the best-fitting results with a lcr scatter given by the dashed lines. Note 
that our model does not include a bulge component. 



sample only the core dispersion, yielding typical values of roughly 
80 km s _1 with a scatter of about 20-30 km s _1 over a range of 
several magnitudes. Unless improved measurements can succeed in 
sampling the disk at larger radii (or at some other means of tracing 
the potential), a proper TF relation may not be easily defined. 

Progress might be made by focusing on the lower, red en- 
velope of the TF relation, rather than attempting to fill in the TF 
plane. The locus of the oldest and reddest points in the TF plane 
generally lie themselves on a well-defined line that serves as the 
high-z/ envelope to the relation. Samples of these old, red objects 
might be easier to construct than complete TF samples, which in- 
clude young, starbursting objects that may have arbitrarily low lu- 
minosities falling off on the younger side of the M peak. SIRTF 
and NGST will have the sensitivity to detect some of the reddest 
objects in the sky, even out to high z. If so, a comparison of the red 
envelope formed by these old objects with predictions, could suf- 
fice to constrain many of the key aspects of galactosynthesis models 
without needing to populate the entire TF plane. Of course, the ef- 
fects of dust will need to be taken into account to so that young, 
dusty systems are distinguished from old galaxies. 

We wish to acknowledge C. C. Steidel and K. A. Adelberger 
for helpful discussions. This work was supported by grants NSF- 
AST-0096023, NSF-AST-9900866, NSF-AST-96 18537, NASA 
NAG5-8506, and DoE DE-FG03-92-ER40701. 



be found on the old side of the M peak. Thus, the TF scatter for 
a fixed mass for such objects would more closely resemble that of 
Figure [l], rather than the more vertically-oriented scatter of Figure 
^| and presumably also lead to a smaller overall TF scatter. The 
differences among the high-z predictions explored here for the TF 
relation and its scatter amount to about 1-2 mag, and 0.2-0.4 mag, 
respectively, so these models should be distinguishable with high 
significance given data of the type in Figure qL The model predic- 
tions for the fjb-M relation can also be used to address the degen- 
eracy. In Figure we see that, owing to the large scatter, the two 
models are difficult to distinguish on the basis of low- z predictions, 
but the tightness of the predicted relation at higher z, if in fact ob- 
served, would make this test a much more powerful discriminator. 



4 DISCUSSION 

We have presented high-z predictions for the TF relation, which 
may help to discriminate among different scenarios for galaxy for- 
mation. The evolution of the spread in the TF relation can probe 
the spread in halo Zf, as well as SF processes in the disk. Specif- 
ically, for observations of evolved systems at redder wavelengths, 
the scatter essentially decouples from the luminosity axis, while 
observations at high z and/or in bluer bands are more sensitive to 
the disk's age, resulting in a scatter that couples more closely to the 
luminosity axis of the TF relation. The predicted [i for a given M 
are fainter and tighter at high z than at low z. 

Our models assume that galaxies in the high-z universe are 
those in the high-z/ tail of the present-day distribution. The circular 
speeds we predict are those for the disk, and so far it has proven dif- 
ficult to extract these reliably from general SF population at z = 1. 
Recent data seem to suggest that beyond the local universe, "nor- 
mal" galaxies actively form stars primarily in their small cores (see, 
e.g., Simard & Pritchet 1998; Rix et al. 1997). Linewidth measure- 
ments of SF galaxies at z = 1 and z — 3 are found to effectively 
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